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Abstract. In this paper, we prove the blowup alternative for Gross-Pitaevskii hier- 
archies on R.^ and give the associated lower bounds on the blowup rate. In particu- 
lar, we prove that any solution of density operators to the focusing Gross-Pitaevskii 
hierarchy blow up in finite time for n > 3 if the energy per some k particles in the 
initial condition is negative. All of these results hold without the assumption of 
factorized conditions for initial values as well as the admissible ones. Our analysis 
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C>J ■ is based on use of a quasi-Banach space of sequences of marginal density matrices. 

(— I , 1. Introduction 

Wh, The time-dependent Gross-Pitaevskii (GP, in short) equation describes the dynamics of 

(-H ' initially trapped Bose-Einstein condensates (see [8, 11] for detailed information). Precisely, 

in units where h = 1 and the mass of the bosons m = 1/2, the condensate wave function 

at time t, (fit = ipt{r), r £K^, satisfies the GP equation 

(1-1) idtift = -Arift +a\ipt\ ifit 

with the normalization J \ipt{r)\^dPr = 1, where a — SivNa is the coupling constant with 
A*' being the number of particles and a the scattering length of the interaction potential. 
This fact was rigorously proved by Erdos, Schlein, and Yau [9, 10, 12, 13], starting from 
a many-body bosonic Schrodinger equation with a short-scale repulsive interaction in the 

^b ' dilute limit. 

^^ , As for the motivation of the present paper, we give a brief description of the main 

ingredients found in [9, 10, 12, 13]. To this end, we assume that the ground state of 
initially trapped Bose gases exhibits a complete condensation (see [16, 17] for detailed 
information). After instantaneously removing the trap, the evolution of A'' bosons is given 
by the Hamiltonian 



JV JV 

(1.2) ffiv = ^(-A,^) + ^l^^(r. 



JH ' The interaction potential is scaled as Viv(r) — N^V{Nr), where 1/ is a positive, spheri- 

. y.' cally symmetric, compactly supported, smooth potential. Let ipN,t be the solution to the 

Schrodinger equation 

(1.3) i9f!/)jv,t = HNiiN,t with ?/'jv,o = ip. 

For fc = 1, . . . , A'^, let Xn t denote the fc-particle marginal density of ipN,t with a consistent 
normalization Tr7jY( = 1. The time evolution of these marginal density matrices is then 
governed by a hierarchy of A'^ coupled equation, commonly known as the Bogoliubov-Born- 
Green-Kirkwood-Yvon (BBGKY, in short) hierarchy 

k k 

(1.4) i9,7^',l = El-^^. ' ^n\] + iN-k)J2 Tr.+i [T^iv(r, - r,+0, 7^'!'^], 
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where Tr^+i denotes the partial trace over r^+i. Furthermore, they show that {iJv t}iv>': 

(k) 

has at least one w -limit point 7^0 t i^i the space of trace class operators, and thereby the 

( k) 

time evolution of 7^0 t i^ described by the Gross-Pitaevskii hierarchy 

k 



(1.5) idtj^Z = Et-^^. ' ^-?'] + ^-s 



loo.t 



J=l 



with the collision operator B^''' — X],=i B, being defined according to 

(1.6) Bf7('=+^'=Trfe+i[5(r,-rfe+0,7*'+"], i = 1, • • • , fc, 

and with cy = J V{r)d'^r. Finally, they proves that for any factorized initial condition 

(1.7) -tt]o = \'Po){'fio\'^\ k>l, with ^oeH\R^), 

the solution of the GP hierarchy (1.5) is unique and factorized: The family 

(1.8) 'y^^]t = \^t){vtf\ fc>i, with ^teH\R'), 

is a solution to (1.5) if and only if (ft is a solution to the GP equation (1.1). 

While the existence of factorized solutions can be easily obtained, the proof of unique- 
ness of solutions of the GP hierarchy is the most difficult part of this analysis and was 
originally obtained in [10] by use of highly sophisticated Feynman graph expansion meth- 
ods inspired by quantum field theory. Their uniqueness result (Theorem 9.1 in [10] or 
Theorem 4 in [11]) reads as: Given a family of density operators F — {'y^''')k>i such that 
(1.9) 
|||7<'='j||. :=Tr[|(l-ArJi•■■(l-A,J^W(l-ArJ^■■•(l-A,J^|]<C^ Vfc > 1, 

for some C > 0, then for any T > 0, there exists at most one solution of density operators 
F(t) = (7<'='(f))fc>i to(1.5) withF(O) = F and such that |||7''''(^)lllfc < C'' holds uniformly 
int e [0,T). 

An alternative method for proving uniqueness has been subsequently proposed by 
Klainerman and Machedon in [15] based on use of space-time bounds on the density 
matrices and introduction of an elegant "board game" argument. Roughly speaking, they 
proved the uniqueness of solutions to (1.5) in the product space "H^ = (^^j^ life of ho- 
mogeneous Sobolev spaces Hj, — H^(R'^'° x R^'°) for all fe > 1, under the assumption of a 
particular a priori space-time bound on the density matrices: 

(1.10) lls^7<'+''LJHl<c^ j = i,...,k, 

with C independent of k (see Section 2 for notions used here). Based on energy con- 
servation, K. Kirkpatrick, B. Schlein, and G. Staffilani in [14] proved recently that the 
inequality (1.10) is indeed satisfied in the case of E'^. 

Let S)^ be the space of all sequences F — (7''°')fe>i of trace class operators such that 



-f{^>0: f:iFlll7'''lll^<l} 



Note that F = (7''°')fe>i satisfies (1.9) if and only if F £ S)^ . Then the uniqueness theorem 
of Erdos, Schlein, and Yau stated above can be reformulated as follows: The initial problem 
for the GP hierarchy (1.5) admits at most one solution of density operators in S^^ . This 
note also applies to the uniqueness theorem of Klainerman and Machedon [15] along with 
the condition (1.10). 

Motivated by this observation, the first named author in [6] introduced a class of 
quasi-Banach spaces for studying the Cauchy problem of the GP hierarchy (1.5). Briefly 
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speaking, a certain Sobolev type space H' of sequences of trace class operators is defined 
as a quasi-Banach space consisting of all F — {'y^''')k>i G ^'tLi ^l such that 



2J Tfe II7 IIh;; < 00 for some A > 0, 



A* 
where H^ = H''(R''" x R*""), equipped with the quasi- norm 



fc=i 

js/jrnkn , , jn)kn\ 



|r||^, — inf<! A> 



0:El^llT''='lk<l| 



It was then proved in [6] that the Cauchy problem for the GP hierarchy (1.5) is locally 
well posed in H" for s > max{i, -^^^j. The precise statement of this that we will use later 
is presented in Section 2. 

In this paper, we continue this line of investigation. First of all, we present some 
notations and preliminaries in Section 2. Then, in Section 3 we will prove the blowup 
alternative for GP hierarchies on R" and give the associated lower bounds on the blowup 
rate. The conservation of energy and Virial identities of solutions to the GP hierarchy 
are presented in Section 4. In Section 5, we prove that any solution of density operators 
to the focusing GP hierarchy blow up in finite time for n > 3 if the energy per some k 
particles in the initial condition is negative. Finally, in Section 6, we will extend all the 
results for the cubic GP hierarchy to the so-called quintic one. 

All these results hold without the assumption of factorized conditions for initial values. 
Yet we do not require the so-called admissible condition proposed in [5] (see Remark 4.2 
below), and hence make some improvements of the associated results found there. 

2. Preliminaries 

As follows, we usually denote by a; = {x^,...,x") a general variable in R" and by 
Xfc = (a;-i, . . . , Xk) a point in R*"" = (R")*. For any x,y £R" we denote by x-y = Yl"=i ^^V^ 
and \x\^ =x-x. For any Xfc,yfc G R*", we set (xfc,yfe) := Yl^^^Xj-yj and |xfcp = {xfc,Xfc). 

We consider the Cauchy problem (initial value problem) for the Gross-Pitacvskii infinite 
linear hierarchy of equations on R", of the form 

r (ift + A('=')7('='(x,;x^)=M[B<''(7r+'')](xfe;x',), 
[7tWo(xfc;Xfe) =% '(xfc;Xfe), fc = 1,2, ..., 

where t e R, Xfe = (xi, X2, ■ ■ ■ , Xk),x.'^. — {x'i,x'2, . . . , x',.) € R*", /i — ±1. Here, 

k k 

A«=^(A.^. -A,,) and b('='=^b('=', 

where Ajj^. refers to the usual Laplace operator with respect to the variables Xj £ R" and 
the operators B^ — B^ ^ — B- _ are defined according to 

[B«(7"=+^')](x.,x',) 

= / dxk+idx'k+i5{xk+i - x'k+i)5(xj - Xk+i)'y'-'''^^\:x.k+i,:x.'k+i), 

and 

= / dxk+idx'k+i5{xk+i - x'k+i)5{x'j ~ Xk+i)-y'-''^^H^k+i,-x.'k+i)- 
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A sequence of functions r(t) = (7^ )fc>i is said to be a Gross-Pitaevskii (GP, in short) 
hierarchy, if they satisfy (2.1) and are symmetric, in the sense that 



7t'°'(xfc,Xfe) =7<'''(x'j^,Xfe) 
and 

(k) / I / \ (k) / f / \ 

7t {xi,...,Xk;xi,...,Xk) =^1 (a;<T(i),---,a;CT(fc);a;CT(i),-- •,a;<^(fc)) 

for any a G Ilfc (life denotes the set of permutations on k elements). 

Let (p e H^(R"), then one can easily verify that a particular GP hierarchy, i.e., a 
solution to (2.1) with factorized initial datum 

k 

7t=o(xfe;x'fc) = Yl<p{xj)ifi(xr), k = l,2,..., 
j=i 
is given by 

fe 
7t'''(xfc;xfc) = Ylv>t{xj)(pt{x'j), fc = 1,2, ..., 
j=i 
where ipt satisfies the cubic non-linear Schrodinger equation 

(2.2) idt^pt ^ —A(pt + ulifit] Vt, (fit^o ^ V^ 

which is defocusing ii jj, = 1, and focusing if /i = — 1. We refer to [1] and references therein 
for the nonlinear Schrodinger equation. 

In the following, unless otherwise specified, we always use ■y^''' , p^''' for denoting sym- 
metric functions in M'=" x K*". For fc > 1 and s e R, we denote by HJ = H''(R'=" x R*^") the 
space of measurable functions 'y^'"'' = 'y^'"\xk,x'^.) in L^(R'°" x R*"), which are symmetric, 
such that 

(2.3) |l7^*^llH= := ||S«7'"||l2(m-xk-) < ~, 

where 

k 

3 = 1 

with the convention S'"'"' — S[ ■ Evidently, H^ is a Hilbert space with the inner product 

/ (fc) (fc)\ / o(k) (fc) rY(fc) (fc)\ 

Moreover, the norm || • ||h» is invariance under the action of e'' , i.e., 

II itAC") (fc)ii II (fc)ii 

||e 7^ '||h| = ||7^ 'I|h| 

because e'' commutates with A^^ for any j. 

Definition 2.1. (cf. [6], Definition 2.1) For s G R we define 

^ oc 00 ^ ^ 

n' = ] (7"=')fc>i G (g)m : E a^II^*''||h| < °«M some\> I , 
I fc=i fc=i J 

equipped with the quast-norrn 

(2.4) j|(7"=')fc>i||^. :=inf|A>0: £ ^fc II^*"IIh| < l| • 

Remark 2.1. Note that || • Wn' is not actually a norm in the sense that it does not satisfy 
||Ar||«s = \^\\\^\\w in general. However, it indeed satisfies the triangle inequality and 
hence, T-i" is a quasi-Banach space. 
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Definition 2.2. (cf. [6], Definition 2.3) For an interval J C R, we define Ll^jH" to be 
the space of all strongly measurable functions r(i) = {7( }fe>i on I with values in H" 
such that 

oc 
ft— i 

equipped with the quasi-norm 



(2.5) 



X^Tfc / II^''°'||h|'^*< oo for some \>0, 

fc=i •'I 

i-norrn 

nmLl^.n^ := inf I A > : ;£ ^ y^ hf ' l|H|rfi < 1 1 



Remark 2.2. For any r(t),$(i) £ Llf^iH" one can verify tliat 

(1) l|Ar(t)||^i^^«. < max{l, |A|}||r(t)|Li^^«. for any A G C; 

(2) ||r(i) + $(t)||^i^^«. < ||r(t)|Li^^^. + ||$(t)||^i^^„.; 

(3) ||r(f)||^i f^s < max{f, |/|}||r(f)||o(7_^3), wliere |/| denotes the Lebesgue measure of 
/; 

(4) l|r(t)|Li ^, = ||r(t)||^i „. + ||r(i)||^i ^. for subintervals h,l2 C / satisfying 

I = hU h and hnh = 0. 



Recall that, in integral formulation, (2.1) can be written as 

(2.6) 7f ' = e"^''\« + f'ds eiC-^JA'-^'^W^Cfc+i)^ ^ ^ ^^ 2, 



/o 

whereafter B^''' — —ifiB^'"'. As noted in [3], such a solution can be obtained by solving 
the following infinite linear hierarchy of integral equations 

(2.7) bW^('=+i) = B('='e"^''+"7('=+^' + f ds B^'^e'^'-'^^'"^"' B^'+'^'ji''+^\ 

Jo 

for any fc > 1. If we write 

'fe>i' 



Ar-(AW7«)^^^ and BT ~ {B<-''^j<-''+^^^ 



then (2.6) and (2.7) can be written as 

(2.8) r(f) = e"^r„ + f ds e'^'-'^^Br{s) 

Jo 

and 

(2.9) BT{t) = Be"^ro + / dsBe'<*"'''^Br(s), 

Jo 

respectively. 

Let us make the notion of solution more precise. 

Definition 2.3. A function r(t) — (7^ )fe>i : / 1— > H" on a non-empty time interval 
G / C R is said to be a local (strong) solution to the Gross-Pitaevskii hierarchy (2.1) if 
it lies in the class C{K, H") for all compact sets K C I and obeys the Duhamel formula 

(2.10) 7f'=e"^"'7i*,''-iM rdse'(*-»)^<'''B('='7l'=+^\ Vi G /, 

Jo 

holds in H^ for every k = 1,2, . . . . 

We refer to the interval I as the lifespan of r(i) with the initial value r(0) = (7g )fc>i. 
We say that r(f) is a maximal-lifespan solution on / = (— Tinin,Tinax) with Tmax = 
rmax(ro) G (0, 00] and Tmin = Tmin(ro) G (0, cx)] if the solution cannot be extended to any 
strictly larger interval. (—Tmin, Tmax) is said to be the maximal lifespan of T{t) with the 
initial value r(0) — (7q )fe>i. We say that r(f) is a global solution if (—Tmin, Tmax) = R, 

I.e., -t niax — -tmin — OO. 
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Remark 2.3. Let •f £ H''(R") and set for A: > 1, 
An immediate computation yields that 

||(7o'°^)fe>l||„= =2||v5||He(K„). 

Thus, for T > Q,ipt e C{{-T, T),W) is a solution to (2.2) with the initial value •ft\t=o = V' 
if and only if 

k 

(2.11) r(t) = (7i")fe>i withif (Xfc;x',) ^\{vt{x,)lM^) 

j=i 

is a solution to (2.1) in C((-r, T), •«=) with r(0) = (7,$'°')fc>i. This yields that the Cauchy 
problem (2.1) in W is equivalent to the one (2.2) in H* for the case that initial conditions 
are factorized. 

The local well-posedness for the GP hierarchy (2.1) in T-i" was proved in [6] for s > 
max{l/2, (n — l)/2}. We state it as Theorems 2.1 and 2.2. 

Theorem 2.1. (cf. [6], Theorem 2.1) Assume that n > 1 and s > ^. The Cauehy problem 
(2.1) is locally well posed. More precisely, there exists a constant A„,s > depending only 
on n and s such that 

(1) For each To = (7o''')fe>i G 'H^ let I = [~T,T] with T = ij^fpr- Then there exists a 

solution r(f) = (7( )k>i G C{I,'H'') to the Gross-Pitaevskii hierarchy (2.1) with the 
initial data To such that 

(2.12) |ir(f)||c(i,Ke)<2||ro||wa. 

(2) Given Jo = [-To, To] with To > 0. // r(t) and r'(i) m C{Io,W) are two solutions to 
(2.1) with initial conditions rt=o ~ To and Fj^q = To *" ^^" respectively, then 

(2.13) ||r(f)-r'(t)||c(/,H») <2||ro-r;,||H3, 

with I = [~T,T], where 

J = mm < If) 



\\r(t)-r'{t)\\cii„.,nn 

Remark 2.4. Theorem 2.1 shows that one has unconditional local wellposedness in "H" 
for the Cauchy problem of the GP hierarchy (2.1) for any s > n/2. This agrees with the 
case of the GP equation (2.2) (see e.g. [19], Proposition 3.8). Using the classical argument 
(see e.g. the proof of Theorem 1.17 in [19]), one can easily shows that given any Fq G W , 
there exists a unique maximal interval of existence I and a unique solution Ft G C{I, T-i"). 
Moreover, if I has a finite endpoint T, i.e., T = Tmax < cxa or T = Tmin < oo, then the 
T-i'-novra of Ft will go to infinity as t ^- T. (Indeed, we will give the corresponding blowup 
rate in Theorem 3.1 below.) Thus, the maximal lifespan of Ft is necessarily open. 

Theorem 2.2. (cf. [6], Theorem 2.2) Assume that n > 2 and s > -^^-S-^. Then, the Cauchy 
problem for the Gross-Pitaevskii hierarchy (2.1) is locally well posed in T-L" . More precisely, 
there exist an absolute constant A > 2 and a constant C = Bn,s > depending only on n 
and s such that 
(1) For every Fo = (7o''')fe>i G H^ let I = [-T,T] with T = „i^"-| . Then there exists 

II'- II ^s 

a solution F(f) = {j'-''\t)) k>i G C{I,'H') to (2.1) with the initial data r(0) = Fq 
satisfying 

(2.14) l|BF(t)]|.i_^, <4A]]Fo]|He. 
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(2) Given Iq = [-To, To] with To > 0. // r(t) G C(Jo,H") so that Br{t) e Ll^i„W 
is a solution to (2.1) with the initial data r(0) = To, then (2.14) holds as well for 
I = [-T,r], where 

r = min<^ To, 



tela / 

(3) Given /o = [-To, To] with To > 0. //r(t) and T'{t) m C{Io,'H'') with Br{t),Br'{t) £ 
Tteio'^" are two solutions to (2.1) with initial conditions r(0) = To and r'(0) = Fg 
in "H^ respectively, then 

(2.15) ||r(f) - r'(f)||c(7.Ha) < (1 + AA)\\ro - r;,||«., 

TOf/i I = [-T,T], where 



T = min < To, 



J^n.s 



B[r{t)-r'{t)]\\l, ^^ + \\ro-r'Ji,. 



In particular, the above results hold for H^ in the case n = 3. 

Remark 2.5. As shown in Theorem 2.2, for the case n/2 > s > {n — l)/2 we require a 
priori assumption Br{t) G LJ^jH" in both the stability and uniqueness parts, although 
we prove that for the existence part, such a priori assumption is not required. At the 
time of this writing, the question remains open whether the condition Br{t) £ LJi^jH" is 
necessary for the uniqueness of solutions. 

Remark 2.6. Note that a similar but different solution space was previously introduced 
by Chen and Pavlovic [3] for studying the initial problem of (2.1). Precisely, given < 
5 < 1, set 

f' oc oo ^ 

w| = r = (7"=').>i e (g)m : lirikj := E^'IIt'^'IIh^ < oo . 

I fc=l fc=l J 

Then, equipped with the norm 11 • ll-^s, "Hi is a Banach space. The local well-posedness 
obtained in [3] states that for any initial data Fq G "H^^ with ^i > 0, there exists a unique 
solution r(i) G C{[-T,T],'Hl^) for some < ^2 < Ci and T > 0, under an additional 
assumption on BT{t) (see also [7] for some improvements in the case s > ^). That is, 
there are two different parameters ^1,^2 in their result. On the other hand, the norm 
||(/Jt||Hs is not compatible with ||r(t)||^s for factorized hierarchies F of the form (2.11). 

As shown in Theorems 2.1 and 2.2, these two undesirable issues will be eliminated if 
the space H" is involved instead of H|. This shows that the space "H" with the quantity 
(2.4) seems more suitable for studying the Cauchy problem of the GP hierarchy (2.1). 

3. Blowup alternative and blowup rate in finite time 

In this section, we prove two results concerning the blowup alternative (e.g., see Def- 
inition 3.1.5 in [1]) and the associated lower bounds on the blow-up rate of solutions to 
the GP hierarchy (2.1). 

Theorem 3.1. Assume that n > 1 and s > ^. If F(t) is a solution to the Gross- 
Pitaevskii hierarchy (2.1) with initial condition r(0) G "H" such that Tmax < 00, then 
limt/iT„ax l|r(i)||-H= = 00, and there exists a constant An^s > depending only on n and s 
such that the following lower bound on the blowup rate holds 

(3.1) l|r(i)||«. > — ^^, vo<i<r„,a>c. 

-t max L/ 

The similar results hold for Tmin- 
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Proof. Suppose Tmax < oo. Wc only need to prove (3.1). For any < t < Tmax, let 
7^t7'*=>(r,Xfc;x'fc) — 7^'=) (^t + r,Xfc; x'^) , r £ R. 

We note that if r(t) = {7''='(i)}fe>i solves the GP hierarchy (2.1), then 7^t^(r) = 
{TZt'y^''\r))k>i is also a solution of the GP hierarchy (2.1). 
To prove (3.1), we fix t £ (0,Tmax). Clearly, 

(3.2) l|7^t7'''(r)||H| = ||7*''(t + r)||H|, 
which implies that 

(3.3) ||7^t^(r)||„e = ||^(i + r)||H=. 
Let r = in (3.3), we obtain 

(3.4) |t7^,^(o)|I«= = ||^(t)||«^<oo. 

Then, Theorem 2.1 implies that there exists r* — An^s/\\'R-tr{0)\\-H= such that 

(3.5) ||7^t^(r)||„. <2||%r(0)||«= 

for any < r < r*. Therefore, combining (3.3) with (3.5) we get 

||r(f + r*)||„. = ||7^t^(r*)||„a < 2||%r(o)||«a = 2||r(t)||«a < oo. 

Hence, the upper maximal-lifespan time Tmax is bounded from below by 

-^ max ^ t I 7" , 

and thus 

A A 

'^^'^ -||%r(o)||«. -||r(i)||„.- 

Consequently, we have 

l|r(t)||T^a > J^"''_ , VO<t<Tn,ax, 
-^ max f- 

as required. 

If Tmin < oo, we cau proceed the same argument and omit the details. □ 

Remark 3.1. Note that in the factorized case, ||r||ws = 2||v5||hs if F = {\ip){ip\'^ )k>i- 
Then the lower bound (3.1) on the blow-up rate coincides with the known one for solutions 
to the GP equation (2.2) in the cases n = 1, 2 ( see e.g. [1]). 

As shown in Theorem 2.2, we require an additional assumption that BF(i) G ijgro t]^" 
for the uniqueness part in the case s < n/2. Then the lifespan / of a solution T{t) should 
be such that 

l|r(i)||c(/,«=) + ||BF(t)|[^i^^^. <^. 

Accordingly, we define the corresponding upper maximal-lifespan time Tmax by 

Tmax = sup{r > : ||F(t)||c([o,T],K=) + II Br(t) || ^1^ ^^_^^ «e < (»}. 

Similarly, we can define Tmin as well. 

Theorem 3.2. Assume that n > 2 and s > -^^-S-i. // F(f) is a solution of the Gross- 
Pitaevskii hierarchy (2.1) with initial condition F(0) G li" such that Tmax < oo, then 
limt/iT„ax l|r(^)ll'H= = OO) o.iT'd there exists a constant Bn,s > depending only on n and 
s such that the following lower hound on the blowup rate holds 

(3.6) l|r(t)||H^ > •^"■° , , VO<i<Tmax. 

\-L max f' j ^ 

In particular, the above inequality holds for the H^-norm in the case n = 3. 
The similar results hold for Tmin. 
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Proof. It suffices to prove (3.6). To this end, we let Tmax < oo and fix T e (0,Tmax) such 
that 



Clearly, (3.2) implies that for any r > 
It follows from (3.4) that 



l|sr(t)|Lj^^^^^^^^„. = l|7^TB^(r)|L.^^^^^^,. 



||7^T^(o)||H. = ||r(r)|iH= < 



CO. 



Then, Theorem 2.2 implies that there exists r* — 77= — ":° ,■> — such that 

(3.7) ||B7^T^(r)|Ll .„. <4A||7^T^(0)|k». 

-re[0,-r*] 

Since 

l|7^T^(r)||c([o,.*].H=) < ||7^T^(0)||7^. + ||B7^T^(r)||^l ^^,, 

it follows from (3.7) that 

(3.8) \\nTT{T)\\cao,r']),n^ < (l + 4A)||7^T^(0)||«.. 
Note that 

\\BT{t)\\ri w» = l|Br(t)|Ui -Hs + ||Br(f)||ri „., 

by (3.7) we get 

(3.9) \\Bmhl^^„^,^^,^n^ - Wmm.l^^^^^^n^ + ll^TBr(r)|L.^_^^^._,. 
By (3.3) again, we obtain 

(3.10) ||r(f)||c([0,T+r.],H») < ||r(t)||c([0,T].-H=) + |17^T^(r) ||c([0,.*],«») 

Also, by the definitions of the operators B and TZt we have 

(3.11) TZTBr{r) = B7^T^(r). 
Then, it follows from (3.7), (3.9) and (3.11) that 

ll^^(*)ll^^e[o,T+.*,«= < ll^r(i)|Li^^^^^«. +4A||7^T^(o)||«. 

(3.12) = l|Br(i)|Li^^^ ^^^. + 4A\\r(T)\\ns 

< 4A(||Br(f)|Li^^^^^^, + ||r(i)||c([o.T].H=)). 

On the other hand, it follows from (3.8) and (3.10) that 

mmc(io,T+r'inn < ||r(t)||c([0.T],«») + (l + 4A)||r(T)||H3 

^ ■ ' <(2 + 4A)||r(i)||c([o,T],«»). 

Thus, from (3.12) and (3.13) we get 



mt)\\cao,T+r'],nn + l|Br(t)||^i^^^_^^^^^^. 



< (2 + 8A)(||r(t)||c([o.T],«») + l|Br(t)||^i^^^_^^„.) < ^. 



^tG[0,T] 

Hence, the upper maxiraal lifespan time Tmax is bounded from below by 

Tmax > T + r* 

and thus 



-L max 1 ^ T — 



|7^T^(o)||?,. l|r(T)|||,. 
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Consequently, we have 

||r(T)||H. > ^'"- 

(r„,ax-r)5 

with Bns = y/Bn,a, ss required. 

The proof for the case of the lower maximal-lifespan time Tmin is similar and omitted. 

D 

Remark 3.2. In the factorized case, the lower bound (3.6) on the blow-up rate coincides 
with the one for solutions to the GP equation (2.2) in the case n = 3 ( see e.g. [1, 2]). 

4. Conservation of energy and Virial identities 

In this section, we will mostly work in Fourier (momentum) space. Following [10], 
we use the convention that variables p,q,r,p' ,q' ,r' always refer to n dimensional Fourier 
variables, while x, x' ,y, y' , z, z' denote the position space variables. With this convention, 
the usual hat indicating the Fourier transform will be omitted. For example, for fe > 1 
the kernel of a bounded operator A on //^(R*") in position space is ii'(xfe; x^), then in the 
momentum space it is given by the Fourier transform 

i^(qfe;q',) = (7^, e-'<--'>'=>e'<-''''^>) ^ J dK,dK',K{^,;^',)e-'<-'-^''^ e'«-<\ 

with the slight abuse of notation of omitting the hat on left hand side. Here, 

k 

(Xfc,qfe) = "^Xj -qj, Vxfc = (a;-i,...,a;fe),qfe = {qi,...,qk) £ R'"" . 
j=i 

Thus, on kernels in the momentum space B'*^' acts according to 

r r>{k) (fc + l)l / ' \ 

= X1 / '^I'^ + 'L'^lk + l 

X ['y'''''^^\pi,---,Pj -qk+1 + q'k+i, ■ ■ ■ ,Pk, qk+i;p'k, q'k+i) 

(4.1) - I'-'^^^Hpk, Qk+i; p'l, ■ ■ ■ , p'j + qk+i - q'k+i,---,p'k,qk+i)j 

= ^ rfqfc+irfqfc+i [ n '^(P' ^ qi)S(pi - qi)\ 

X 7*'°^^'(qfc+i;q'fc+i){'5(Pj -q'j)5{pj - fe +17^+1 -q'k+i]) 

- Hpj - q3)S{p'j - [q'j + q'k+i - gfc+i]) j- 

Following [5], we introduce 

k 

(4.2) Ekir) := X: / dx,dx^5(x, - x^) (Iv.^. ■ V.. 7^" + jB^':W'+'^ 

Note that 

/ dxdx'5{x — x')Wx ■ W'xA{x,x') 

(4.3) = - / dxdx'5{x - x')AxA{x, x') 
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which was previously proved in [5] . Then 

k 

(4.4) Sfe(r) :=^ /'dx,dxU(xfe-x^)(-iA.^.7('=> + Jb]'_?7"=+'' 

3 = 1-^ 

Further, we note that each of the terms in the sum equals to the one obtained for j = 1, 
by the symmetry of 7''°'(xfe; x^) and j^'''^^' {xk+i;x'i._^i). Indeed, since 

v., •V,.7<'='(xfe;x'0 = ^^^j^^7*'''(xfc;xL), 



we have that 



dxfcdxfe5(xfe - Xfe)V^j • V^/7'*'(xfe;x'fc) = j dx^ ^ --^—- j-7^'='(xfe; x^) 



Then, by the symmetry of 'y^^' {pi.k\'^'k) with respect to the components of Xfc and xj. one 
has 

dxfcdxfe(5(xfc - Xfe)V:c3 ■ '^^'7*''' = / dxfcdxfe(5(xfe - Xfe)Va;, • Vj;'7^*^ 



for all 1 < J,i < fc. We note that the calculation for the interaction term was presented in 
[5]. Thus, 

(4.5) £,(r) = fc|dx,dx',5(x,-x',)(iv., •V^7^'' + ^<i7"=+'' 

Remark 4.1. For factorized states F = {\ip){ip\® )k>i with <p e H^(E") one finds that 



"''"-'' ''^\\y^\\l^ + ^v-"' 



Edr) ^ kMir'' i^Wv^lli^ + ^y^. 



In this case, Ei(T) is the usual expression of the conserved energy for solutions of (2.2). 

In what follows, we prove energy conservation of k particles for solutions F(t) to the 
GP hierarchy (2.1) for any k. 

Theorem 4.1. Assume that F(t) is a solution of the Gross-Pitaevskii hierarchy (2.1) with 
initial condition F(0) G T-L" . Then Ek(T(t)) is a conserved quantity, i.e., 

(4.6) Ek{V{t)) = Ek{V{Q)), ViG(-r„,in,T^,,). 

Remark 4.2. This result was previously proved in [5] under the assumption that V{t) is 
admissible in the sense that 



(fe) 
Tt 



(xfc;x'fc) = / dxk+i-yt \^k,Xk+i;:x!k,^k+i), Vfc > 1. 



In that case, it can be shown that EkiVit)) = kEi{V{t)) and thus reduces to show that 
Ei{V{t)) is a conserved quantity. On the other hand, we note that for a sequence of "pure 
states" 7'''-'(xfe;x^) = ?/>(xfe)i/)(x5.), the admissibility condition implies that (7'*-')fe>i must 
be factorized, by Schmidt's decomposition theorem (see e.g. [18]). This indicates that the 
admissibility requirement seems a little restrictive. 

Proof. First of all, by (4.5) we have 

Ek(J{t)) = kj d^kd^k5{^k-0(^V., ■ V,,7f' + f Bl'i7r^ 

Then, by (2.1) we have 

idtEk{T(t)) = k[(I) + (II) + (III) + (IV)] 



where 
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(//) : = I Adxfcdx',5(x, - x',)V., ■ v.. B^^^f +^\ 



(III) : = -^ y dxfedx',5(xfc - ^',)B[':iA^'+'hr'\ 
{IV) : = i|dx.dxU(x. ~x',)<iB('=+i)7r+''- 

In order to prove dtEk{r{t)) = 0, we shall prove that (/) = {IV) = and (//) + (///) = 0. 
For term (/), we note that 

-2(7) = / dxfedx'fe5(xfe - x'fe) v., • V,; A('=)7f ' 



(fc) 



= ^ '[{S{XI - x'i)dXkdx'kS{Xj - x'j){Aa:^ - A^/)V:„i • '^x{'yl 

= ^ / ]^5(x! - a;J)dxfcdXfe5(a;-j - x^) 

3 = 1 iT^J 

X /e'«^-P'=>-<^-P'^»dprfp'(lp,f - ip;i')(pi •p'i)7*''(p;p')- 



Since 



j dxidx',dp,dp',5{x, -x',)e'^^y''^-^'r^'^\\p,\' - lp;i')(pi ■p'i)7<''(P;p') 
= / dpjdp'jS{pj -p'j){\pA'^ - |p^l^)(pi ■pi)7*'°'(p;p') = 0, 



for any j, we proves (J) — 0. 
For {IV), we have 

4(/l/) =|dx,dxU(x, -x^)B(;=JB('=+^)7(''+^' 
fc+i 



= Y1 j dxkdx'k5{xk ~ x'k)B['^l(^j['''^^\xk+i,Xj;x'k+i,Xj) 

(fe+2)/ / / /n\ 

-7t (xfe+i,a;j;Xfe+i,Xj)j 

= 5Z / rfxfcdx'fcc5(xfc -x'fc)(^7t**+^'(xfc,a;i,a;j;Xfe,a;i,a;j) 
j=i"^ 

_ (fc+2)/ /, / /n\ 

Tt t,Xfc, a^i, a^j , X;., Xi, Xj- M 
y~^ / dxfc l7j*'°+^'(xfc,a;i,a;j;Xfe,a;i,%) - 'y^''^'^\xk,xi,Xj;Xk,xi,Xj)\ 

.7 = 1 "^ 



J = 

=0 



This proves {IV) = 0. 
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For (//), we have 
(//) =f I dxferfx^5(xfc - x^)V., ■ v.. B('')7r+^' 



=f Y. / d^^d^'.si^, - xL) v., ■ v., [Bl';l - B^';1W+'^ 



k 

2 ■ 

+ f E / d^^d^Ji^k - xL) v., • v.; [B^X^l - B(;!^]7f+i) 

=1 /'dxfedx^5(xfe -x^)V., ■ V.,[<j -<!]7r+'' + 1^ /"dxfedx',5(xfe -x',) 

X [(A^i7t'°"^^')(xfc,a;j;Xfe,a;j) - (A^i7(*'°+^')(xfc, a;j;Xfe,a;^)j . 
Note that 

/ dxfcdxfec5(xfc -Xfe)|^(A:,i7*'°+^')(xfc,a;j;Xfe,a;3) - (A:,i7t*'°+^')(xfe, a;j;Xfe,a;^-)J =0. 
We get 
(//) =1 1 dxfedx^5(xfe - x^)V., ■ v., (B+,+1 - i3i-,+i)7f +'^ 

= - ^ / dxfcdx'fc5(xfe -x'fc)|^(A^/^7(*''+^')(xfe,a;i;x'fc,a;i) - (A^i7j'*+^')(xfe,a;i; x'^, a;i)J 
where we have used the identity (4.3). 



Now we turn to the term [III)- By symmetry of 7} ' in xj; and xj., we have 
(777) = - J fd^kd^,5{^k - x',)73(:=iA('=+i)7.*'=+i) 



j=i 



= - I E / '^Xferf'^'fc'5(xfe - x',)73(;=j(A.^ - A,, )7f +" 



■ / dxfcdx'fc(5(xfc - x'fc)B(*j(A^i - A^. )7j'* 



fc+i 
4' 



i=2 



^^ / dx,dxU(x, -xL)i3^;=j(A.^ - A,,)7f+'' 



I y"dx,dx',5(x, - x',)73« (A,, - A,, )7('=+i' 



= - 2 / dyikdyikS(yik - x'^.) 

X [(A:ci7t'°^^')(^fe'^i''^'fc'^'i) ^ (A^'^7t'°"'"^')(xfe,a;i;x'fc,a;i) 

Thus, (77) + (777) = 0. 

In summary, we have dtE{I'{i)) = 0. Therefore, E{I'{t)) = 75(r(0)) is a conserved 
quantity. □ 

Now we turn to Virial type identities for the GP hierarchy (2.1). This is necessary for 
the application of Glassey's argument for blowup in finite time of solutions to the focusing 
(m = -1) GP hierarchy (2.1). 
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Given a solution r(t) to (2.1), we define for any fc > 1, 

fc 
T4(r(f)):=Tr[|xfe|V"W] =X^ /dXfe|:c,|V'Ht,Xfc;Xfe). 

By tlie symmetry of 7'*^', we have 

\4(r(f)) = fc /'dxfc|xi|'7<'='(f,Xfe;Xfc). 

Then, the following Virial type identity for the GP hierarchy (2.1) holds. 

Theorem 4.2. Assume that r(i) = (7^ )k>i solves the Gross-Pttaevskii hierarchy (2.1). 
Then for any fc > 1, 



(4.7) dfVk{r{t)) = 8k f dpk\pif-(^''\pk;Pk) + 2nk^i f dxk-f^''+'\^k.. 



Xi;Xk,Xi) 



Remark 4.3. The identity (4.7) was previously proved in [5] under the assumption that 
r(i) is admissible (see Remark 4.2 for detailed information). 

Proof. We write 

7''^(x,;x'fe) = |dpfedp',e'<^'-'''=>~'<''^'''^^''>(Pfc;p'fe) 
and define 

p(xfe) — 7^'''(xfe;xfc) = / dpfcdpfee'^P'="'''=''''=^7**'(pfc;p'fc). 
Then, 

a,p(x,)=ydpfedpU'^'"="^''=''"=^a7'"(Pfe;p'0 
(4.8) =iydpfcdpU'^''^"'''^'"''^(A('='7*")(Pfc;Pfc) 

-Mi|dp.rfPfee'<P'=-P''='^'=>(BW7<'+'')(Pfe;p'0- 
We first have 



iydpfcdpU'<P'=-P^--'=>A('='7('='(pfe;p',) 



= -i / dp.rfpU'^^^^^-^'-^dp.l' - |p'.l')7*''(Pfe;Pfe) 
(4.9) ^ 

= -i / dpfcdpfee'^P'' ^''"'''"^(pfc + p'fc,Pfc - Pfc>7*'°'(pfc;pt 



k) 



= -Vx, • y dp,dpU'<^'="^-"'=>(Pfe +p'fc)7'"(Pfe;p'fc). 
On the other hand, for j = 1,2, . . . ,k, 
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= / dqfc+idqfe+ic5(q'fc - p'k)S{qk+i + q^ - q'k+i - Pj)'[{S{pi - gi)7*'°^^'(qfe+i; Qfe+i) 

= / dijfe+i(igfc+i7*'°+^'(pi,.. .,Pj-i,gfc+i - gfc+i +Pj,Pj+i,...,Pfc,<jfe+i;Pfe,'7fc+i) 
Similarly, one has 

B},_7 HP'=;Pfc) 

= |dxfedx',e-'<P'"='^>+'<P'^--''=>B(;=27('=+^'(xfe;x',) 

= |dXfedx',dq,+idqUie-'<'''='"'=>+'<P''='^''=>e'(<'''='^'=>+*+i-^^<'>^'^^>-'^+i-^^'7''+''(qfc+i;qUi) 

k 

dqfc+idqfe+i5(qfe - Pfc)(5(gfc+i - g^ - g^+i +P3)]^5(p| - q'ih''^^^H<ik+uq'^+i) 



1^3 



■- / dqk+idq'k+i-y'-'^^^\pk, qk+i; p'l, ■ ■ .,p'j-i,qk+i - q'k+i +Pj,p'j+i, ■ ■ • ,p'fc,9fc+i)- 



Therefore, 



rfpfedpU'<''^"'''^'"'=>(Bf7<'+'')(Pfe;p'.) 



= I dpkdpie'^^--^'-'-"^ [(B«7<'=+'')(p.;Pfe) - (B]'27<'+'')(Pfe;p'.)] 

= / dgfc+idgfe+idpfcrfpfee'^^'°"^'"'''''^*'°"'"^'(Pi. • • ■,<l'k+i - Qk+i +Pj, ■ ■ ■ ,Pk,qk+i;p'k,<l'k+i) 

- / (igfc+id(7fc+idpfcdpfee'<P'="P'=''''=^<'=+^'(pfc,(7fc+i;pi,...,gfc+i - q'k+i + p'j, ■ ■ ■ ,p'k,<l'k+i) 

= 0. 

where the last equality is obtained by applying the change of variables pj — >■ q'^._^_i—qk+i+Pj 
and p'j —7- q'k+i — qk+i +p'j in the second term of the second equality so that the difference 
Pj — p'j remains unchanged. Since 

we conclude that 

(4.10) |dpfedp',e'<'"=-'''^'"'=^(B<''7^'+")(Pfe;p'.) = 0. 

Therefore, by combining (4. 8), (4. 9) and (4.10) we have 

(4.11) Stp + Vx, -P^O, 
where 



p -.^ j dpkdp'j^^'^-^'-'-'^^pk + p'kW'Hpk^Pk) 



Z. Chen, C. He, and C. Liu 

Now, we define 



The time derivative is given by 

(4.12) dtM = f dxfc(xfe, dtP) = Im + IIm, 

where 

Im = 1 (ixfc(ipfedp'fee'^'''="'''=''''=^(xfc,pfe + p'fe)(A<'°'7<*')(pfe; p'j.), 

and 

IIm = -/ii /"dxfedpfedp^e'<P'=-P'fc-^'=>(xfe,pfe +pL)(B("'7*'+")(Pfc;p'fc)- 
For the term Im we have 
Im = ~i / dxfcdpfedp'fce'^P'="'''=''''=^(xfe,pfc +Pfc>(|pfc|^ - |p'fc|^)7*'''(Pfc; p'fc) 
= -i / dxfcdpfedp'fce'<P'="P'=''''=>(xfc,pfc +Pfe)(pfc +Pfe,Pfc - p'fc>7*'''(Pfc; p'fc) 
= -i / dxfc(ipfedp'fce'*P'="P'=''''=^((pfe +p'fc)(pfc +p'fc)^Xfc,pfe - p'fc>7*'°'(pfc; p'fc) 
= - I rfpfedp',7<"(Pfc; p'k) I rfxfe((pfe + p^)(pfc + p^)^Xfe, Vx,e'<P'=-P'^-^'=>) 
= I dp,dpW\pk;p'k)Tr{p, + p^)(p, + p^)^ I dxfee'<P^-P''='''^> 
= / dpfcdpfe7*'''(pfc;pfc)|pfc +pfe|^5(pfe -p'fc) 
= 4 / rfpfc|pfc|^7*'°'(Pfc;p'fc), 



where x^, Pfe, p'^ are all considered as nfc x 1 matrices, and A'^ denotes the transpose of a 
matrix A. 

Next, we determine the term IIm- To this end, for each j = 1, 2, . . . , A; we have 

rfpfedpU'<''^-''-''^>(p,+pL)(Bf7<'=+'')(Pfe;p',) 

= |dp.dpU'<''^-'''--'=>(p. +pL)[(B]'_?7"=+'')(p.;Pfe) - {B\':li^'+'^){p,;p',] 
dpfcdp'fcd(7fc+id(jfe+ie'^P'="P'"'''=^(pfc +Pfe) 



/' 



X ['-f'''''^^\pi,---,i'k+i ~ Qk+i +Pj, . . . ,pk,qk+i;p'k,q'k+i) 
-7 [Pk,qk+i;pi,. ■■,qk+i - qk+i +Pj,---,Pk,qk+i)\ 

In the last term, we apply the change of variables pj — >■ pj — qk+i + q'k+i s^nd p'j — >■ 
p'j — qk+i + q'k+i so that the difference Pj —p'j remains unchanged. Then the above integral 
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equals 



dpkdp'^:dqk+ldq'^:+le^^''^^>'•'^''''{pk + p'k) 

X i^'^^^Hpi, ■ ■ ■,q'k+i - Qk+i +Pj, ■ ■ ■,Pk,qk+i;p'k,q'k+i) 

- / (ipfedp'fcd(7fc+id(7fc+ie'<P'="P'=''''=>(pi +p[,...,pj + p'^ - 2qk+i +2q'k+i, . . . ,Pk + p'k) 

X 7*'°^^'(pi, • • ■,q'k+i - qk+i +Pj, ■ ■ ■,Pk,qk+i;p'k,q'k+i) 

= / dpkdp'kdqk+idq'k+ie'-^'''^>"'''''>2{0, ..., qk+i - q'k+i,- • • , 0) 

X i^'^^^Hpi, ■ ■ ■,q'k+i - qk+i +Pj, ■ ■ ■ ,pk,qk+i;p'k,q'k+i) 



The contribution of this term to the integral J dxfc{xfc, dtP) is given by 

-^i / dxkdpkdp'kdqk+idq'k+ie'-^''~^>"''''hxj ■ (qk+i - q'k+i) 

X 7*'°+^'(pi,...,(7fc+i -qk+i +Pj, ■ ■ ■ ,Pk,qk+l■,p'k,q'k+l)■ 
Now, we apply Fourier transform again. This integral equals 

(4.13) 
-^i / dy:kdpkdp'kdqk+idq'k+idyk+idy'k+ie'-^'''^>"''''hxj ■ {qk+i - q'k+i) 

X g-i<Pfe.yfc>-i(9fc + i-gfe+l)-!/j-i9fc + l-!/fc+l+i<pL.yL>+i'?fc + r!'fc + l^('=+l)(-yj,^j.yJ, ) 

= ~^li I dxfcdyfc+idyfc+i7'''+^^(yfc+i;yfc+i) / dpkdp'kdqk+idq'k+i 

= -^ / dxfedyfc+idyfc+i7^''+^'(yfc+i;y'fc+i)(5(xfc - yfc)(5(xfe ~ y'k) 
X f dqk+idq'k+i2xj . Vy^.e''"»+i''''^"'"»+i'"'«^+i-(^^-"^^+i^ 

= -P d^kdyk+idyk+i'y^''^^\yk+i;xk,y'k+i)S{yik -yk) 

X 2xj ■ Vy.S{yj - yk+i)5{yj - y'^+i) 
= -^ / dxfedyfe+idyfe+i7*''+^'(xfc,jyfc+i;Xfc,yfc+i) 

X 2a::j -V^^Sixj - yk+i)S{yk+i ~ y'k+i) 
= -^ / dxkdyk+i'y^''^^\xk,yk+i;xk,yk+i)2xj ■'Vxj5{xj - yk+i) 

= ^ / dxkdyk+i5{xj - yk+i){2n + 2xj ■ Vxj)j'-''^^\xk,yk+i;xk,yk+i)- 
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Note that 



dxkXj ■ Va;j7*'°+^'(xfc,a;-j;xfe,a;j) 
= / dTikdyk+i5{xj -yk+i){xj -V^. + yk+i ■ Vy^^j7*'°+^'(xfc, y^+i; x^, y^+i) 

^' ' = dxkdyk+i5{xj -yk+i)[xj ■Vx^'y'-'^^^\xk,yk+i;xk,yk+i) 

+ Vk+i ■ Va^^i7^''+^'(a;i, . . . , yk+i, . . . ,Xk,Xj;xi, . . . , yt+i, ...,Xk,Xj)] 

= / dxkdyk+i5{xj -yk+i)2xj ■ V^^'y'-''^'^'' {xk,yk+i;Xk,yk+i) 

where we used the symmetry of 7'''°^^-', and renamed the variables in the last term. Also, 

(4.15) / dxkXj ■\7xjy'-'^^^''{xk,Xj;Xk,Xj) = -n / dxk'y'-''^^\xk,Xj;Xk,Xj) 

from integrating by parts. 

Therefore, combing (4. 13), (4. 14) and (4.15) yields 

IIm =^A' / dxkdyk+i5{xj - yk+i){2n + 2xj ■'Vxj)'y \xk,yk+\;xk,yk+i) 
j=i ■^ 

= ^/i / dxkdyk+iS{xj ~ yk+i){2n~ n)j^''^^\xk,yk+\;Xk,yk+\) 

=nkfi / dxfe7*'°+^'(xfc,a;i;xfc,a;i). 
Finally, we combine (4.11) and (4.12) to conclude that 

d't / dxfc 1x^1^7**' (f,Xfe;Xfe) 

= 2 f dxk{xk,dtP) 

= 8 / dpfc|pfe|^7*'°'(pfc;pfc) + 2nfc^ / dxfc7*'°+^'(xfc, xi;Xfe,a;i) 



= 8fc / dpk\pifj'-''\pk;Pk) + 2nkfi / dxk'y'-'''^^\xk,xi;Xk,xi). 
This completes the proof of Theorem 4.2. D 

5. Blowup of solutions to the focusing GP hierarchy 

In this section, using conservation of energy and Virial type identities obtained in the 
previous section, we prove a result on blowup of solutions to focusing GP hierarchies in 
finite time. 

Recall that f)^ is the space of all sequences F = (7 ) of trace class operators satisfying 
||F||j5i < 00, where 

l|r|Ui=mf|A>0: ^J^|||7('=)|||, <l|, 

and 1 1 17'*' I life = Tr[|5'('='7('='|] for any fc > 1. 

Theorem 5.1. Let n > 3. Assume that F(f) — (7^ )fe>i solves the focusing (i.e. fi = —1) 
Gross- Pitaevskii hierarchy (2.1) with initial condition F(0) = (79 )fe>i € S)^ . If {ji )fe>i 
is a sequence of density operators such that '14(F(0)) < 00 and Ek{T{0)) < for some 
k > 1, then the solution F(f) blows up in finite time with respect to Sj^ . 
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Remark 5.1. We note that in the factorized case, Theorem 5.1 coincides with the corre- 
sponding result for the focusing GP equation (2.2) (e.g., see Theorem 6.5.4 in [1]). 

Proof. By (4.7), we have 

9?14(r(t)) = 16fcTr[ - iA,,7W] +2nkfiTr[B[';lj<-'^+^^] 

= 16Ek{r{t)) + fik{2n - 4)Tr[B$"J7('=+^>] 

= 16£;fc(r(0)) - k{2n - 4)Tr[B(^J7('=+'>] 

where we have used Theorem 4.1. Since {-yl )k>i is a sequence of nonnegative operators 
and n > 3, we have 

(5.1) 9tVfe(r(t)) < i6£fc(r(o)). 

Since Vfe(r(t)) is nonnegative, we conclude from the assumption Ek(r{0)) < that there 
exists a finite time T* such that Vk{T{t)) \ as t /" T*. 
On the other hand, we have 

l = Tr[7"=)(t)]<(Tr[|xi|%«(t)])'(Tr[^7"='(t)])' 

where we have first used the Cauchy-Schwarz inequality, and then the Hardy inequality 
(noticing that n > 2). Thus, 

Tr[-A.,7'''(t)]=Tr[(-A.j5(_A,,J^W(i)]>_i_^^ as t ^ T* , 

where we have used (4.3) in the first equality. 
Now, by (2.3) and (2.4), one has 



l|r(t)||^,i > |||7<''(t)llk=Tr[5('='7<'='(t)] =Tr[n(l-A.^)7<'='(t) 

j=i 

k 

> ^ Tr [(- A.^. )7<'=' (t)] = fcTr [ - A., 7*'' (i) 
j=i 
as t /^ T* , which establishes blowup in finite time. D 

6. The quintic Gross-Pitaevskii hierarchy 

In this section, we consider the so-called quintic Gross-Pitaevskii hierarchy and present 
the corresponding results similar to the ones obtained for the cubic case in previous sec- 
tions. 

6.1. Preliminaries. Recall that the quintic Gross-Pitaevskii (or, GP in short) hierarchy 
r(t) = (7('='(f))fc>i is given by 

(6.1) ia.7.'''' = [J2{-A..^),jn +^Q''hr'\ M = ±l, 

3=1 

in n dimensions, for fc G N, where the operator Q'*' is defined by 

3 = 1 

It is defocusing if /i = 1, and focusing if /i = —1. We note that the quintic GP hierarchy 
accounts for 3-body interactions between the Bose particles (see [4] and references therein 
for details). 
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In terms of kernel functions, the Cauchy problem for the quintic GP hierarchy (6.1) 
can be written as follows 

l7S(xfc;x',)=7^'='(xfe;x',), fc G N, 

where QC^' := E,tj Qf with the action of Qf> = Q</=j - Qfl on 'y^''+^\^k+2,^'k+2) G 
^(]jj(fc+2)n ^ ]^(fc+2)„^| ^gjjjg defined according to 

(Q«7''=+^')(x.,xL) 

:= / da;fe+ida;fe+2rfa;'fe_|_ida;'fe_|_27*'°"^^'(xfc,a;fc+i,a;fe+2;Xfe,a;fe+i,a;'fe_|_2) 

fe+2 

X TT 5{xj — xi)5{xj — x'l) 

_ (fc+2)/ . / \ 

and 

(Q<:=!7"=+'')(x.,x',) 

:= / dxk+idxk+2dx'^,_f_idx'k+2'y^'"^'^H^k,Xk+i,Xk+2;:>i'k,x'k+i,x'k+2) 

fc+2 

X TT (5(a;j — xi)5{x'j — x'l) 

l=k + \ 
(fc+2) / / / , / / / \ 

I \^k ; Xj , Xj , X^ , Xj , Xj f , 

for j = l,...,k. 

Let (j9 £ H^(R"), then one can easily verify that a particular solution to (6.2) with 
initial conditions 

fc 

7t=o(xfe;x'fc) = Yl<p{xj)<p{xr), fe = 1,2, ..., 
J=l 

is given by 

k 

7f'(xfc;xfe) = ]^v't(a;j)v5t(a;^0> fc = 1,2, ..., 
j=i 
where tpt satisfies the quintic non-linear Schrodinger equation 

(6.3) idttpt — —A(ft + n\ft\ ft, ipt^o — f- 
The GP hierarchy (6.2) can be written in the integral form 

(6.4) ^W^e"^<"7('='+ rdse'(*-^'^<'"0«7f+'', fe = l,2,..., 

Jo 

where Q ~ — i/iQ^'^ . Evidently, such a solution can be obtained by solving the following 
infinity linear hierarchy of integral equations 

(6.5) 0('=>7r+'^ = 0('=)e«^"°+''7('=+^> + /* ds Q^'^U'^'-^^^'"^"' Q^''+^^^i'=+^\ 

Jo 

for any fc > 1 . If we write 

AF := (A('=>7('=>)fc>i and QF := (Q<'='7''+'')fc>i, 
then (6.4) and (6.5) can be written as 

(6.6) F(t) = e"^Fo + f ds e'('-"^^(5F(s), 

Jo 
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and 

(6.7) Qr{t) = Qe"^ro + [ ds Qe'^'-'^^Qr{s), 

Jo 

respectively. 

Formally we can expand the solution 7^ of (6.4) for any m > 2 as 

Tt =e To + Z^ / "*2 / dt4- ■ ■ dt2je ^' Q' ' ■ ■ ■ 

j^i Jo Jo Jo 

X pi(t20-l)-«23)A<'= + "0-l»^(fc+2(j-l)) itjA<'=+23) (fe+2j) 

(6.8) V f^ To 

+ /' dt, r dU--- /''"""'' dt2™e''*-*^'^'''Q(''' • • • 
Jo Jo Jo 

with the convention to = t- 

6.2. Local well-posedness. In this subsection, we present the local well-posedness re- 
sults for the quintic GP hierarchy (6.2). 

Let us make the notion of solution more precise. 

Definition 6.1. A function V{t) — (7^ )fe>i : I 1— > T-L" on a non-empty time interval 
G / C R is said to be a local [strong) solution to the Gross-Pitaevskii hierarchy (6.2) if 
it lies in the class C{K, H") for all compact sets K C I and obeys the Duhamel formula 

(6.9) 7.*"' = e"^*''7('=' - i/i /' ds e'^'-'^^"" Q^'^ji'+'\ Vt £ /, 

Jo 

holds in HJ for every k = 1,2, . . . . 

The following theorem was stated in [6]. 

Theorem 6.1. (cf [6], Theorem 5.1) Assume that n > 1 and s > ^. The Cauchy problem 
(6.2) is locally well posed. More precisely, there exists a constant K„^s > depending only 
on n and s such that 

(1) For every To = (7,'i*')fe>i & W , we let T = ..^"^J and I = [-T, T]. Then there exists 

a solution r(t) = (74 )fe>i £ C{I,'H") to the Gross-Pitaevskii hierarchy (6.2) with 
the initial data To satisfying 

(6.10) ||r(t)||c(/,«»)<2||ro||H. 

(2) Given lo = [-To, To] with To > 0. //r(f),r'(f) m C{Io,'H'') are two solutions to the 
Gross-Pitaevskii hierarchy (6.2) with the initial conditions r(0) — To and r'(0) = Fq 
inW, respectively, then 

(6.11) ||r(t) - r'(t)||c.(/,«.) < 2||ro - rj,||„. 

with I = [-T,T], where 



^---^^-M^)-r'(.)|i^,,,,., 

The proof can be obtained by use of the fully expanded iterated Duhamel series and a 
Cauchy convergence criterion, based on the following inequality 

||qW^('=+2) 11^. < c„,.fc||7<'=+'' ||h|^, , Vfc > 1, 

with Cn,3 > being a constant depending only on n and s, which was proved in [4] 
(Theorem 4.3 there). 

For the case s < n/2 we have 
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Theorem 6.2. Assume that n > 2 and s > ^^^- Then, the Cauchy problem for the 
Gross-Pitaevskii hierarchy (6.2) is locally well posed in T-L" . More precisely, there exist an 
absolute constant A > 2 and a constant C — M„^s > depending only on n and s such 
that 

(1) For every To = (7,'j''')fe>i G W, we let T = *^"J and I = [-T,T]. Then there exists 

a solution r(t) = (7*'''(i))fc>i G C{I,'H") to (6.2) with the initial data r(0) = To such 
that 

(6.12) IIQr(t)|Li^^„. <2A||ro||«.. 

(2) Git;en Jo = [-To, To] with To > 0. If r{t) G C(/o,'H") so i/ioi Qr(f) G Ll^i^W is 
a solution to (6.2) wii/i the initial data r(0) = To, then (6.12) holds true as well for 
I = [-T,T], where 

rj. ■ [^ M^.s 1 

i = mm \ Jo, 



IIQr(t)r „. + ||ro|||,. 



(3) Git;en /o = [-To, To] u-ii/i To > 0. lfV(t) and V {t) m C{Io,'W') with QT{t),QT'{t) G 
Ll^i^'H" are two solutions to (6.2) with initial conditions r(0) — To and r'(0) = V[) 
inW, respectively, then 

(6.13) ||r(i) - r'{t)\\c(i,nn < (1 + 2^)liro - Folk-, 

with I = [-r,T], where 

T = mm < To, 



\\Q[rit)-rit)]\\l ^. + ||ro-r;,r«. 



This result was announced in [6]. For the sake of completeness, we present the proof. 
To this end, we need two preliminary results 

Lemma 6.1. Assume that n > 2 and s > ^^-^. Then there exists a constant Cn,s > 
depending only on n and s such that, for any symmetric ^'-^+^'^ £ 5(R**+^'" x R'''+^'"), 

fr 1 A\ II^C^) itA<'= + ^' (fe+2)|| ^ r-i II (fc+Z) I 

(6-14) \\Q) 'e 7^ + lli?(H^) ^ C„,s\W ^ '\Wi^^ 

for all k>l, where j — 1,2, . . . ,k. 

Consequently, Q^^' can be extended to the space HJ^2 such that 

(6.15) IIO<'''e'-<^"\('=+^)|L.,H|) < C„,.fc||7"=+^'||H|,, 

/ora;;7('=+2)GH^+2. 

This lemma was proved first for n = 3 in [15], and then in [3] (Proposition A.l there) 
for general case. 

For any F — {pji )k>i we define 

P.+2„(F)(i):= fdt, rdt,--- r""'dt2.e'('-'^'^"'+^'Q('=+^)... 
Jo Jo Jo 

with the convention t = to- 

The following lemma is crucial for the proof of Theorem 6.2. 

Lemma 6.2. Assume that n > 2 and s > ^^^. Then there exists an absolute constant 
A > 2 and a constant Cn,a depending only on n and s so that the estimates below hold 



Gross-Pitacvskii hierarchies 



(1) For any To = (7i*,")fe>i G 0^=1 H^ 

(6.16) ||Q"='Pfe+2..(e"^ro)(t)Lj^^H| < fc^'+^(c„..r)^|i7r''+''llH=^,^+, 

for k,j > 1 andT >0, where I = [-T, T] and e^^Fo = (e"^"'''7o''')fe>i- 

(2) For anyT>0 and r(i) = (7t''°')fe>i mi/i ji''^ £ Ll^i_j,.j,]'Hl, 

<fcA'=+™(C„,sT)^||(5'*=+'''"'7(''+'""+2^(i)||ii ^H| , , 
for k,m>l, where I = [-T,T]. 

Proof. The inequalities (6.16) and (6.17) can be proved by using the so-called "board 
game" argument presented in [15]. For the details see the proof of Theorem 6.2 in [4] 
(e.g., Proposition A. 2 in [3]). D 

Now we are ready to prove Theorem 6.2. To this end, we introduce the system 

(6.18) r(i) = e^'^ro + [ ds e'('""'^H3 







(6.19) Et = Qe"^r„ + / ds Qe'C-^'^H 



which is formally equivalent to the system (6.6) and (6.7). 
The proof is divided into three parts as follows. 

Proof (1). Let n > 2 and s > {n ~ l)/2. Let To = {7o''''}fc>i G W and Ho = 
{po (^)}fc>i — 0. Given A: > 1, for any m > 1 we define 

(6.20) p«(i) = Q"='e"^""^''7('=+^' + r ds Q«e'('-^'^"=^''pL'+?'(-), 

Jo 

for t G / = [-T, r], where T > will be fixed later. Set Em{t) = {pm'(t)}fc>i for every 
m > 1. By iterating (6.20), for every m > 2 one has 

m — 1 

pL"(t)=E^*''^^+2.:'(e"^ro)(t) 

with the convenience Pfe+2,o(e'*^ro)(i) ~ e''^ 7o • By Lemma 6.2 (1) we have 



E<3*"-P^+2..(e"^ro)(i)|| ^ 



j—0 *£^ fc 

m — 1 



I.e., 



(6.21) IIpL'^IL;,,h| < k Y: A'+^{^cZTy^'\hr''^''\\nu 



j=0 
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Set T ■- „^'";{i with M„,, = w^-^- For A > one has by (6.21) 

'-Oils oOtis/I -j \ / 



'A 

fc>i 



E-'- II (fc)ii 

. 771 — 1 

fe>l j=0 

fe>l £=fe+2 



< 



fe>i 

'A\\VoU^\k _ 1 e 



^E'- 



A^^ V A / ^(1-0^' 



fe>l 



with ^ = A||ro||H»/A. Note that A > 2, choosing ^ = i, in other words, A = 2A||ro||T.«= 
we have 

El II ik)\\ ^ 1 i ^ -I 

fc>l ^ ^' 

This concludes that for every m, > 1, Hm G L\^iT-C and 
(6.22) ||H™||^i^^^e<2A||ro||K.. 

Now, for fixed fc > 1 and any n, m, with n > m we have 

11 (fc) (fe)ii 

, n-1 



<^^^.+,(y^;^).+i||^(.+2.+2)| 



j=m 



f>fc + 2m+4 ^" " ■' 

This concludes that for each fc > 1, pm converges in L^gjH^ as m — )> oo, whose limitation 
is denoted by p^ \ 

Set H(f) = {p<'''(i)}fc>i- Note that for any m > 1, 

II rrf.QWe'(*-^)-<^"^'[pL^+?)(.)-p('=+^)(.)lj| ^ 

<E r r rfirf^iiQ^.^e'(-«)-'^"^'[pL^+?'(.) -p('=+^)(.)]ii^. 

<T^/^^ r d.l|Q,,e'(-»)-<^"^'[pL^+f (.) -p('=+^)(.)]|L. „. 



<r i-T^/^ii (fc+2) .(fe+2)|| 

where we have used the Cauchy-Schwarz inequality with respect to the integral in t in 
the second inequaUty and used Lemma 6.1 in the last inequality. Thus, taking m — > oo 
in (6.20) we prove that H is a solution to (6.19). Moreover, taking m — > cx) in (6.22) we 
obtain (6.12). 
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(2). Fix To > 0. Suppose r(f) £ C{Io,W) is a solution to (6.2) so that QV{t) e 
^telo T \W ■ Let T G (0, To], which will be fixed later. By iterating (6.5), for every m > 2 
one has 

m — 1 

+ Y. 0*"^fe+2.i(e"^ro)(t) + Q("'Pfe+2,„(r)(i). 

By Lemma 6.2 we have 






^J+l|U,('=+2i+2)| 

j=0 

I 7 tk-\-m / fT^ ;T;\m||^(fc + 2m) (fc + 2m + 2)|| 

Set 

J = mm < Jo 



Qr(i)||*, ^, + l|ro| 



4 



Taking m — )• oo we have 






i^llTo IIhs 



Then, for A > we have 

1 
A 






< 



^^V A ) Z^(||ro|k.)^ll^» IIh|^ 



fc>i ^=1 ^" """ ^ ^ ^' 

with ^ = A||ro||«s/A. Note that A>1, choosing A = 2^||ro||ws, that is, ^ = i, we have 

^A^ll^ 7* iLi H= ^ ^(i_^)2 ^^- 

This completes the proof of (2). 

(3). Fix To > with Jo = [-To, To]. Suppose V(t)Xii) S C(Io,W) are two solutions 
to (6.2) such that QV(t),QV'(t) e L\^i„U' with r(0) = Po and r'(0) = Fq in W, 
respectively. Since (6.2) is linear, it suffices to consider F(f) instead of T{t) — r'(t). Set 

^ . ]rr M^,s 1 

T = mm < To, 



IIQr(f)r «. + l|ro|||,. 

Then, by (2) we have for / = [-T, T], 

l|r(t)||c(i,H.) < IIFoIIk. + |lQF(i)||^i^^^. < (1 + 2A)||Fo|lwe. 
This completes the proof. □ 

6.3. Blowup alternative and blowup rate in finite time. By slightly repeating the 
proofs of Theorems 3.1 and 3.2 with the help of Theorems 6.1 and 6.2, respectively, we 
can obtain two results concerning the blowup alternative of solutions to the quintic GP 
hierarchy (6.2) as follows. 
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Theorem 6.3. Assume that n > 1 and s > ^. IfT{t) is a solution of the Gross- Pitaevskii 
hierarchy (6.1) with initial condition r(0) G W such that Tmax < oo (resp. Tmin < oo), 
then linit/'Tniax l|r(i)||-Hs = oo (resp. liint\,-T„i„ ||r(i)||-M» ~ oo^, and the following lower 
bound on the blowup rate holds 

l|r(i)||«e > — ^"'' 1 , VO<i<rn,ax, 

(6.24) ^ \Tm..-t\2 

|T„.in-t|t 

where Kn,s > is a constant depending only on n and s. 



0|r(-t)||«a > ^^^^-r, VO<t<Tn,in), 



Theorem 6.4. Assume that n > 2 and s > ^^^-^. // r(i) is a solution of the Gross- 
Pitaevskii hierarchy (6.1) with initial condition r(0) G ?^^ sttc/i i/ioi Tmax < oo (resp. 
Tmin < 00^, then limt^Tmax l|r(i)||Ks ~ oo ('resp. limt\,_T„i„ ||r(i)||K= = oo^, and the 
following lower bound on the blowup rate holds 

l|r(i)||K= > ^"-^ 1 , VO<t<Tn,ax, 

(6.25) iTma^il' 

(l|r(-t)||T,a > ^^"-^ , , VO<t<Tn,in), 

J niin (-1 * 

where Mn,s > is a constant depending only on n and s. 
We omit the details of the proofs. 

6.4. Blowup in finite time. In this subsection, we present energy conservation and 
Virial type identities of solutions to (6.2). As an application, we can obtain a result on 
blowup in finite time of solutions to the focusing GP hierarchy (6.2). The arguments 
are the same as those involved in Sections 4 and 5. The required modifications are not 
difficult and left to the interested reader. However, for the sake of convenience, we write 
the corresponding definitions and results. 

First of all, as in [5] we introduce the energy functional of solutions to (6.2) as follows. 

E,im) := iTr[^(-A.^.)7"='] + ^Tr [X: Q^:=i7"='] 

3=1 j=l 

for any fc > 1. Then 

Theorem 6.5. Assume that r(f) is a solution of the Gross- Pitaevskii hierarchy (6.2) with 
initial condition r(0) G W . Then Ef^{T(t)) is a conserved quantity, i.e., 

(6.26) Ek{r{t)) = Ek{r{0)), ViG (-Tmin, Tmax). 
To state the Virial type identity for solutions to (6.2), we set 

14(r(i)) - Tr[|xfe|S'''] = J2 /rfxfclx.f 7"=)(i,Xfe;Xfe). 

Then we have 

Theorem 6.6. Assume that r{t) — {-yl )k>i solves the Gross-Pitaevskii hierarchy (6.2). 
Then for any k > 1, 

(6.27) dtVk{r{t)) = Sk f dpfe|pi|S''='(pfc; pfe) + |nfc^ f dxfc7('=+''(xfe, a;i; x^, xi). 

By Glassey's argument, we can conclude from Theorems 6.5 and 6.6 the following 
result. 
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Theorem 6.7. Let n > 3. Assume that r(f) — (7J )fe>i solves the focusing (i.e. fi = —1) 
Gross-Pitaevskii hierarchy (6.2) with initial condition r(0) = (7,) )fe>i £ f)"'. If {'ft )fe>i 
is a sequence of density operators such that \4(r(0)) < 00 and _Efe(r(0)) < for some 
fe > 1, then the solution r(t) blows up in finite time with respect to Sj^ . 
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